arXiv: 1504.00725v2 [cond-mat.mes-hall] 4Jun2015 


Spin splitting in 2D monochalcogenide semiconductors 


Dat T. Do, 1 ’* Subhendra D. Mahanti , 1 and Chili Wei Lai 1 ’^ 

1 Department of Physics and Astronomy, Michigan State University, East Lansing, MI f8824, USA 

(Dated: June 5, 2015) 

We report ah initio calculations of the spin splitting of the uppermost valence band (UVB) 
and the lowermost conduction band (LCB) in bulk and atomically thin GaS, GaSe, GaTe, and 
InSe. These layered monochalcogenides appear in four major polytypes depending on the stacking 
order, except for the monoclinic GaTe. Bulk and few-layer t- and 7 -type, and odd-number /3-type 
GaS, GaSe, and InSe crystals are noncentrosymmetric. The spin splittings of the UVB and the 
LCB near the T-point in the Brillouin zone are finite, but still smaller than those in a zinc-blende 
semiconductor such as GaAs. On the other hand, the spin splitting is zero in centrosymmetric bulk 
and even-number few-layer /3-type GaS, GaSe, and InSe, owing to the constraint of spatial inversion 
symmetry. By contrast, GaTe exhibits zero spin splitting because it is centrosymmetric down to 
a single layer. In these monochalcogenide semiconductors, the separation of the non-degenerate 
conduction and valence bands from adjacent bands results in the suppression of Elliot-Yafet spin 
relaxation mechanism. Therefore, the electron- and hole-spin relaxation times in these systems with 
zero or minimal spin splittings are expected to exceed those in GaAs when the D’yakonov-Perel’ 
spin relaxation mechanism is also suppressed. 


I. INTRODUCTION 

Potential applications in spin-dependent electronics 
and optoelectronics have driven the search for materials 
capable of exhibiting a high degree of spin polarization 
and long spin relaxation time 1 ’ 2 . However, optical gener¬ 
ation of electron and hole spin polarization and resulting 
polarized luminescence are typically limited by the mix¬ 
ing of degenerate valence bands in most semiconductors 2 . 
Recent reports of valley polarization in atomically thin 
transition metal dichalcogenides (TMDs ) 3-7 suggest po¬ 
tential exploitation of both spin and valley degrees of 
freedom for electronics and optoelectronics. In an exper¬ 
imental study 8 , we demonstrated the high generation and 
preservation of optical spin polarization and dynamics in 
a group-III monochalcogenide, GaSe, under nonresonant 
optical pumping. The observed near unity optical spin 
polarization 9,10 is attributed to suppressed electron and 
hole spin relaxation rates resulting from reduced valence- 
band mixing. However, the microscopic spin relaxation 
mechanisms in GaSe and related monochalcogenides are 
not fully understood. 

In metals and semiconductors, the major spin relax¬ 
ation mechanisms- including Elliott-Yafet (EY ) 11,12 and 
D’yakonov-Perel’ (DP ) 13-15 mechanisms- are associated 
with the spin-orbit interaction (SOI) and the spin-orbit- 
induced spin splitting, A s (k) = \E(k,^) — E{k,\f ) | 2,16 . 
Considering spin-relaxation with a four-state (two bands 
with spin) model Hamiltonian in the absence of an ex¬ 
ternal magnetic field, one can relate the spin relaxation 
rate of electrons (holes) when the Fermi energy (corre¬ 
sponding to Fermi vector ) is away from the conduc¬ 
tion (valence) band edge with the following equation 16 : 
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where T p = h/r p is the scattering rate of the elec¬ 


tron/hole, with Tp being the corresponding momentum 
scattering (or correlation) time, A s (fc) being the spin- 
orbit-induced spin splitting, and L(k) being the SOI be¬ 
tween the adjacent bands with energy separation A g . In 
GaSe, the p 2 -like uppermost valence band (UVB) is well 
isolated from the lowermost conduction band (LCB) (~2 
eV) and the adjacent p^-like valence bands, and as a re¬ 
sult L/A g r; 0.02-0.04 17 ’ 18 . The hole-spin relaxation due 
to the EY mechanism Tf ' 1 ss (-^^Tp is thus expected 
to be much smaller than the momentum relaxation rate 
Tp. The spin relaxation caused by the DP mechanism 
can be seen as being due to the precession of spins in 
an effective magnetic field associated with A s (fc) 2 ’ 15,16 . 
The DP spin relaxation rate is proportional to the spin 
splitting, oc TpA s (fc) 2 , where t p is the momentum 

relaxation time. Therefore, when the spin relaxation is 
dominated by the DP mechanism, the smaller the spin 
splitting, the longer the spin relaxation time r s = h/T s 
for the same momentum relaxation rate T p . 

To understand the spin relaxation, one first needs the 
momentum(/c)-dependent A s (k) of the bands near the 
fundamental gap. In the absence of magnetic fields, 
A s (k) is zero in centrosymmetric crystals because of the 
constraints of time-reversal symmetry [E{k, j') = E(—k, j, 
), Kramers degeneracy] and spatial inversion symme¬ 
try [E{k, t) = E(—k, j')]. When the inversion sym¬ 
metry is broken in crystals (bulk inversion asymmetry 
(BIA )) 19 or heterostructures (structural inversion asym¬ 
metry (SIA)) 20 ’ 21 , A s (k) is finite, and only the Kramers 
degeneracy is left. Understanding A s (k) in GaAs and 
other zinc-blende semiconductors has been a subject 
of considerable interest 22-27 since the seminal work of 
Dresselhaus 19 . Ab initio calculations, such as LDA (or 
GGA) and self-consistent GW methods, of A s (k) in 
bulk GaAs and two-dimensional GaAs-based superlat¬ 
tices and heterostructures have improved the understand¬ 
ing of the spin splitting 23-26 . A few theoretical calcula- 
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tions of A s (k) in TMDs also have been reported 28,29 . In 
this study, we report ab initio calculations of A s (k) of 
the uppermost valence band (UVB) and the lowermost 
conduction band (LCB) in GaSe and related group-III 
monochalcogenides, including GaS, GaTe, and InSe. 



(b) 




FIG. 1. (a) Side view of the 2Ha e-, 3R 7 -, 2Hb 0-, and 4H 
(5-polytype MX (M = Ga, In; X = S, Se) unit cell, (b) Top 
view of the MX single layer, (c) Side view of the monoclinic 
GaTe unit cell, (d) Top view of GaTe single layer. M and 
X are big (brown) and small (green) spheres, respectively. In 
the centrosymmetric systems, one possible inversion center is 
denoted by a red circle. 

Monochalcogenides MX (M = Ga, In; X = S, Se) 
crystallize in hexagonal layered structures 30 (Fig. 1) of 
four major polytypes, namely e, 7 , 0, and S (Fig. la), 
depending on the stacking order (hereinafter referred to 
as MX crystals), e-, 7 -, 0-, and 5-MX crystals be¬ 
long to the space group (Schoenflies notation) of D\ h , 
Cf,,, Dq H , and Cq v , respectively. Monolayer MX crystals 
(space group D^h) are noncentrosymmetric. Bulk e-, 7 -, 
and 5-MX crystals, which appear in an AB, ABC, and 
ABCD stacking order, are noncentrosymmetric, while 0- 
MX crystals are centrosymmetric with an AB stacking 
order. Additionally, there are two exceptions: (1) an 
atomically thin 0-MX crystal with even-number layers is 
centrosymmetric, and (2) a bilayer 5-MX crystal can be 
identical to either a bilayer e-MX crystal (noncentrosym¬ 
metric) or a bilayer 0-MX crystal (centrosymmetric) de¬ 
pending on which two layers are isolated from a bulk 5- 
MX crystal. GaTe appears as a distorted form of the 
MX structure, where one out of three Ga-Ga bonds lies 
in the a-b plane (Fig. lc-d). In contrast to MX crystals, 
GaTe crystals belong to the monoclinic lattice system 
(space group (7^), and are centrosymmetric down to a 
single layer 30 . 

Bulk and few-layer e- and 7 -type, as well as odd- 
number few-layer /3-type GaS, GaSe, and InSe crystals 
finite spin splittings, while bulk and even-number few- 


layer /3-type GaS, GaSe, and InSe as well as GaTe crys¬ 
tals exhibit zero spin splitting. The difference is due to 
the constraints of the aforementioned time-reversal and 
spatial inversion symmetry (or the lack of it). 


II. COMPUTATIONAL METHODS 

We compute the band structures and A s (k) of va¬ 
lence and conduction bands with the projector aug¬ 
mented wave method as implemented in the VASP 31 33 
package and the full-potential (linearized) augmented 
plane-wave as implemented in the WIEN2k 36,37 package. 
The band structures are calculated with the WIEN2k 
package, with the optimized crystal structures deter¬ 
mined by minimizing the total energy with all electrons 
(including core electrons) with VASP. In all calcula¬ 
tions, exchange-correlation energies are determined by 
the Perdew-Burke-Ernzerhof (PBE ) 38 generalized gradi¬ 
ent approximation (GGA) 39 , which systematically un¬ 
derestimates the band gaps and produces E(k) disper¬ 
sions (effective masses) different from experimental val¬ 
ues. These shortcomings of the GGA also limit the ac¬ 
curacy of the calculated A s (fc). 

The spin-orbit interaction (SOI) is included in our cal¬ 
culations of the overall band structure and the spin split¬ 
ting of a given band in a self-consistent manner using 
a second variation approach 40-42 . The SOI Hamiltonian 
in the spherical symmetric potential can be represented 
as: H so = 1/(2m 2 c 2 ) x [1/r dV(r)/dr\L ■ 5, where m e 
is the electron mass, c the speed of light, L and S the 
orbital and spin momentum vectors, and V (r) an effec¬ 
tive single particle local potential seen by the electron. 
This form of H so is correct as long as V(r ) is local and 
isotropic. In Hartree approximation and LDA, the ef¬ 
fective potential is indeed local, though it is not always 
isotropic. The isotropic approximation is valid because 
the dominant contribution to H so is from regions near 
the nucleus. However, local approximations do not give 
correct band structure near the band gap. The accuracy 
of the band gap can be improved with hybrid models such 
as the Heyd-Scuseria-Ernzerhof (HSE06) 43 (a mixture of 
non-local and local exchange) or GW-like theories. In 
these approximations, the aforementioned simple form 
of H so is not necessarily valid 44,45 . A compromised ap¬ 
proach invloves calculating the band structure using non¬ 
local or GW theories while calculating the spin splitting 
with a local form of the H so and the LDA potential 23,46 . 
Nonetheless, one has to examine uncertainties in these 
hybrid calculations critically. 

To model a few-layer thin film, we create a supercell 
(supercell method 39 ) containing one few-layer structure 
and a 15-25 A thick vacuum spacer, which is large enough 
to suppress interactions arising from the artificial peri¬ 
odicity present in the supercell method. The crystalline 
c-axis of the supercell is set perpendicular to the crys¬ 
talline a-b plane. In this way, one can distinguish the 
effects of intra- and inter-layer interactions on the elec- 
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tronic structures in few-layer structures. The number of 
atoms in a unit cell is as follows: eight for e- and /3- 
MX, twelve for 7 -MX, sixteen for 5-MX , and twelve 
for monoclinic GaTe. To obtain an energy accuracy of 
0.1 meV in self-consistent calculations, we use T-centered 
Monkhorst-Pack 47 fc-meshes of 24 x 24 x 4 and 24 x 24 x 1 
for bulk and few-layer GaSe-type structures, respectively. 
For GaTe, we use meshes of 16 x 6 x 8 and 16 x 6 x 2 for 
bulk and few-layer GaTe, respectively. 


III. RESULTS 

A. Band structure: bulk versus single-layer 



FIG. 2. Electronic band structures, along K-r-M in the 
hexagonal Brillouin zone (BZ), of bulk (solid blue curves) and 
monolayer (dotted red curves) (a) /LGaS, (b) e-GaSe, and (c) 
e-InSe. The zero energy is set at the valence band maximum. 

In Fig. 2, we show the electronic band structures of 
bulk and monolayer /3-GaS, e-GaSe, and e-InSe, which 
are the most naturally abundant. The general features 
of the electronic band structures, except the spin split¬ 
ting, are nearly polytype-independent, owing to the weak 
inter-layer interactions. The lowermost conduction band 
(LCB) has s-like symmetry, whereas the two uppermost 
valence bands (UVBs) have p z -like symmetry. The p x ,y~ 
like valence bands appear ~1 eV below the UVB as a 
result of the crystal field and SOI. The calculated band 
structures for e-GaSe show a nearly direct band gap at 
the T-point of the Brillouin zone (BZ), where a valley ap¬ 
pears in the UVB. The energy of the LCB at the T-point 
is ^0.5 eV lower than that at the M point, consistent with 
the hybrid density functional calculations 48 . On the con¬ 
trary, tight binding calculations show that the energy of 


the LCB at the M point for GaSe is ~10 meV below that 
at the T-point in the BZ 49 . 

The band gap is seen to decrease with increasing 
atomic number (Ga—Un or S—>Se). The calculated band 
gaps are 2.0 eV, 1.3 eV, and 0.71 eV for /3-GaS, e-GaSe, 
and e-InSe, respectively, which are each smaller than the 
experimental values (~3.1 eV, 2.0 eV, and 1.3 eV) 30 . The 
band-gap underestimation can be remedied with, for ex¬ 
ample, the HSE06 hybrid functional 43,48,50 . In the ab¬ 
sence of SOI, the p x>y states are doubly degenerate at the 
T-point. On the other hand, the SOI lifts this energy de¬ 
generacy with a spin-orbit splitting A so ~ 0.09 eV, 0.34 
eV, and 0.31 eV in GaS, GaSe, and InSe, respectively. 
A so in GaSe and InSe are similar in magnitude, but a 
factor of three smaller in GaS, agreeing with previously 
reported calculations 51-54 . A so in GaS is minimal, as ex¬ 
pected from the weak SOI in the lighter S anions which 
govern the characteristics of the few uppermost valence 
bands of GaS. Monolayer GaS, GaSe, and InSe have very 
similar band structures (Fig. 2). We note two different 
features in the band structures of monolayer MX s in 
comparison with their bulk counterparts: ( 1 ) the quan¬ 
tum confinement along the c-axis increases the band gap 
to 2.36 eV, 1.78 eV, and 1.4 eV for GaS, GaSe and InSe, 
respectively, and ( 2 ) the band gap becomes indirect as 
the valley at the T-point becomes wider in momentum 
(k = |£|) and deeper in energy (E). 


Rrillrmin yrme 



FIG. 3. The electronic band structure of monolayer GaTe 
(left) along the selected high-symmetry directions in the 2D 
BZ (right). 

The band structure of GaTe (bulk) has also been cal¬ 
culated with GGA 55,56 , showing a direct band gap of ~1 
eV. The inclusion of SOI causes negligible changes in the 
UVB and LCB of GaTe. Monolayer GaTe shows a direct 
band gap of 1.4 eV (Fig. 3), with LCB having two nearly 
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degenerate minima at the T and C points. At the C point, 
LCB has s-like symmetry while the UVB has p y -like sym¬ 
metry. SOI removes the p x , y degeneracy of valence bands 
at r, with A so ~ 0.2 eV. A so is smaller in GaTe than 
in GaSe despite Te being heavier than Se. The reduc¬ 
tion in the strength of A so is due to the quenching of 
orbital angular momentum in the lower symmetry crys¬ 
talline structure, as demonstrated by a sizable A so ~ 
0.7 eV calculated for a hypothetical /3-type GaTe (space 
group D\ h ). 


B. Spin splitting 



than 1 meV, and LCB spin splitting is smaller by a factor 
of five compared to that of the monolayer. The thickness 
dependent spin splitting in /3-GaSe presented here are 
consistent with those reported in MoS 2 29 , which has the 
same symmetry as /3-GaSe. Bulk y-GaSe has similar spin 
splittings as bulk e-GaSe, with decreasing spin splittings 
as the number of layers increase. 



- - - GaAs (GGA) - ■ - GaAs (GW) 


-GaSe GaS InSe 

FIG. 5. Spin splitting A s (k') along the T-K line for (a) the 
uppermost valence band (UVB) and (b) lowermost conduc¬ 
tion band (LCB) of monolayer GaSe, GaS, InSe, and bulk 
GaAs. For valence bands of GaAs, we show the splitting for 
the heavy hole (HH), which is the UVB, the light hole (LH) 
and the split-off (SO) bands. The spin splittings calculated 
with the GW method are extracted from Refs. 23,26 . 


FIG. 4. Spin splitting A s (k') as a function of k' = k/kmax 
(kmax = TK) for the uppermost valence band (UVB) and the 
lowermost conduction band (LCB) along the T K in n-layer 
e-GaSe (a-b) and j3- GaSe (c-d) (n=l,2,3, and oo (bulk)). 

In Fig. 4, we show the spin splittings of the UVB 
(A v s (k)) and the LCB (A c s (k)) along the T-K direction 
in e- and /3-GaSe. The spin splitting along the T-M di¬ 
rection is zero, obeying the constraint of spatial inver¬ 
sion symmetry. Both A v s {k') and A c s (k') decrease with 
the number of layers, approaching those in the bulk. At 
k' = k/km ax = 0.15 (kmax is k at the K point in the 
BZ), Ag ss 6 meV and 4 meV for monolayer and bulk 
e-GaSe, respectively. The nearly layer-independent LCB 
spin splitting has a value A^ ss 11 meV at k' = 0.15, 
which is slightly larger than A^. 

In contrast to e-GaSe, bulk and even-number few- 
layer /3-GaSe crystals have zero spin splitting (Fig. 4c- 
d), obeying the constraint of spatial inversion symmetry. 
The Ag(/c') and A c s {k') in odd-number few-layer /3-GaSe 
crystals are finite, but diminish rapidly with increasing 
layers. In trilayer /3-GaSe, the UVB spin splitting is less 


In Fig. 5, we compare A v s (k') and A c s {k') in monolayer 
GaS, GaSe, and InSe (group-III monochalcogenides) and 
bulk GaAs (a representative zinc-blende III-V semicon¬ 
ductor). Among the monolayer group-III monochalco¬ 
genides, overall spin splittings decrease from GaSe, to 
InSe, and then to GaS. The spin splittings typically in¬ 
crease with the increasing atomic number of constituent 
atoms as result of the enhanced SOI in the heavier atoms. 
However, other details of the band structure such as the 
band gap also contribute to the spin splittings. 

The valence band of GaAs consists of a heavy hole 
(HH), a light hole (LH), and a split-off (SO) band 22,27 . 
The calculated HH spin splitting is close to that in the 
UVB of GaSe. However, the spin splittings in the LH 
and SO bands are at least a factor of two larger than 
that in the UVB of GaSe. The calculated overall LCB 
spin splitting in GaAs is also larger than that in GaSe. 
The magnitude at k' = 0.15 is more than two times larger 
in GaAs than in GaSe. The spin splitting of the heavy- 
hole band is reduced by about a factor of two when the 
GW method is used in lieu of the GGA method. 
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IV. DISCUSSION 


Band 

UVB 

LCB 

Coefficient 

A (meV) 

B (eV) 

A (meV) 

B (eV) 

Monolayer 

1.0 

4.9 

0.3 

4.7 

2 -layer 

2.9 

2.6 

2.0 

4.2 

3-layer 

0.3 

3.0 

0.4 

4.8 

4-layer 

0.6 

2.6 

0.4 

4.9 

5-layer 

0.3 

2.5 

0.4 

4.9 

6 -layer 

0.4 

2.4 

0.4 

4.9 

Bulk 

0.1 

2.2 

0.2 

5.0 


TABLE I. Linear (A) and cubic ( B ) coefficients of the In¬ 
dependence of spin splitting, A „(k') = Ak ' + Bk’ 3 with k' = 
k/k max < 0.05. 


The spin splittings discussed above concern mainly the 
overall spin splitting up to k' = 0.15. To understand the 
spin relaxation mechanisms, we need to identify the k- 
dependence of the spin splitting in the vicinity of the 
T-point. At small k , the k ■ p theory predicts that, 
in noncentrosymmetric zinc-blende and wurtzite struc¬ 
tures, the fc-dependence of the spin splitting contains 
both a linear and a cubic term when the core levels are 
considered 19,26,57,58 . To illustrate the ^-dependence of 
the spin splitting, we fit the calculated A s (k') in e-GaSe 
with the function A s (k') = A k' + B k' 3 for k' < 0.05 
(Table I). The energy scales for the coefficients A (meV) 
and B (eV) are consistent with those determined from 
GW calculations for GaAs 23,26 . Although B is three to 
four orders of magnitude large r than A, there exists a 
crossover value of k' = i JA/B ~ 10~ 2 below which the 
linear term dominates. In contrast to the GaAs case 
where the linear term is negligible for the LCB, we find 
a sizable linear term for the LCB in GaSe. The cubic 
coefficient B for the LCB is ~4-5 eV for the monolayer 
to the bulk, with the bilayer case being slightly differ¬ 
ent. In contrast to the UVB, there appears to be an 
odd-even-layer effect: B values for the odd layers (1 and 
3) are larger, but are close to the bulk values for the even 
layers (2 and 4). For the LCB, the A values are similar 
for all the layers, except for the bilayer (A = 2.0 meV) 
and the bulk (A = 0.2 meV). Note that bilayer GaSe 
has an unusually large A value, and for the UVB, there 
appears to be an odd-even effect like that in the LCB. 
The A value for the bilayer is nearly three times that 
for the monolayer, whereas A for the four-layer is two 
times that for the trilayer. As pointed out in the case 
of GaAs, these subtle differences are due to the charac¬ 
teristics of the UVB and LCB energy values and wave 
functions, and their mixing with other bands including 
the core levels 19 ’ 26,57 . 

The DFT-based theories such as LDA (or GGA) un¬ 


derestimate band gaps and do not give accurate effective 
masses, resulting in overestimated A s (fc) 23,25 ’ 26 . GW 
calculations reproduce more accurate band parameters, 
such as the band gap and effective mass, but are compu¬ 
tationally more intensive than LDA (GGA) calculations. 
For simplicity, in this work, we have used GGA to cal¬ 
culate A s (k). The GGA calculation underestimates the 
GaAs band gap by a factor of ten; however, the spin 
splitting only deviates from that determined by the GW 
calculation by a factor of two. The band gaps are under¬ 
estimated by the GGA for GaSe and related monochalco- 
genides by a factor of approximately two, which is signifi¬ 
cantly less than that for GaAs. For example, in GaSe, the 
GGA calculation gives a band gap of about 1 eV, which is 
off from the GW/HSE06 calculation 50 and the measured 
band gap (~2 eV) 51 by a factor of two. Therefore, we ex¬ 
pect the GGA calculation to produce spin splittings close 
to the value obtained with the GW calculation. We also 
expect similar variations of A s with k from one conduc¬ 
tion/valence band to another and from bulk to atomically 
thin layers. 

V. CONCLUSION 

We present a systematic study of spin-orbit-induced 
spin splittings bulk and atomically thin group-III 
monochalcogenides MX' (M = Ga, In; X' = S, Se, 
Te). The spin splitting vary with anion element and 
crystal symmetry. Centrosymmetric crystals, including 
bulk /3-type GaS, GaSe, and InSe, as well as mono¬ 
clinic GaTe down to the monolayer, have zero spin split¬ 
ting, as anticipated from the constraints of spatial in¬ 
version symmetry and time-reversal symmetry. In these 
monochalcogenide semiconductors, the separation of the 
non-degenerate conduction and valence bands from other 
adjacent bands results in suppression of Elliot-Yafet spin 
relaxation mechanism. Therefore, the electron and hole 
spin relaxation times in these systems with zero or mini¬ 
mal spin splittings and reduced valence-band mixing are 
expected to be longer than those in a zinc-blende semi¬ 
conductor (eg., GaAs 22,27,59 ), owing to the suppression of 
D’yakonov-Perel’ and Elliot-Yafet spin relaxation mech¬ 
anisms. 
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